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ABSTRACT 

The internal rotation of post-main sequence stars is investigated, in response to the convective 
pumping of angular momentum toward the stellar core, combined with a tight magnetic coupling 
between core and envelope. The spin evolution is calculated using model stars of initial mass 1, 
1.5 and 5 M e , taking into account mass loss on the giant branches. We also include the deposition 
of orbital angular momentum from a sub-stellar companion, as influenced by tidal drag along with 
the excitation of orbital eccentricity by a fluctuating gravitational quadrupole moment. A range of 
angular velocity profiles Q(r) is considered in the envelope, extending from solid rotation to constant 
specific angular momentum. We focus on the backreaction of the Coriolis force, and the threshold 
for dynamo action in the inner envelope. Quantitative agreement with measurements of core rotation 
in subgiants and post-He core flash stars by Kepler is obtained with a two-layer angular velocity 
profile: uniform specific angular momentum where the Coriolis parameter Co = flT con < 1 (here T con 
is the convective time); and D(r) oc r -1 where Co > 1. The inner profile is interpreted in terms 
of a balance between the Coriolis force and angular pressure gradients driven by radially extended 
convective plumes. Inward angular momentum pumping reduces the surface rotation of subgiants, 
and the need for a rejuvenated magnetic wind torque. The co-evolution of internal magnetic fields 
and rotation is considered in Kissin & Thompson, along with the breaking of the rotational coupling 
between core and envelope due to heavy mass loss. 

Subject headings: stars: giants - stars: white dwarfs - stars: rotation - magnetic fields 


1. INTRODUCTION 

Measurements of internal oscillation modes of red giant 
stars by the Kepler satellite have been used to infer rapid 
rotation in the radiative core ■ in some cases, a n or der of 
magnitude faster t han at the surface ()Beck et al.ll2012t 
iMosser etT al. 2012) . Thi s evidence has been interpreted 
by iCantiello et all (12014H in terms of a decoupling of the 
rotation of the core from the convective envelope, which 
is assumed to rotate as a solid body. 

Here we reconsider the role that convection plays in 
redistributing angular momentum within a star. So far 
only a limited understanding has been developed of the 
interaction between rotation and convection. The most 
convincing theoretica l calculations involve d irect numer¬ 
ical simulations (e.g. iBrun fc Palaciosl l 2009 ). 

Red giants offer a promising test case, because (i) the 
influence of the Coriolis force is reduced (although not 
entirely eliminated) in comparison with main sequence 
(MS) stars; (ii) the great depth of the convective envelope 
reduces the effect of boundary conditions on the angular 
velocity profile; and (iii) feedback from magnetic fields 
on the rotation profile is less important, especially near 
the tip of the red giant and asymptotic giant branches 
(RGB and AGB), where most of the envelope rotates 
slowly. Disentangling the most important processes that 
control the rotation profile of a star is, for these reasons, 
easier in red giants than in the Sun. 

The internal rotation of giant stars also has important 
implications for the magnetism and spin of their white 
dwarf (WD) remnants. A basic theme of this paper, and 


the companion paper (iKissin fc Thompson! 12015L here¬ 
after Paper II) is that one must consider the co-evolution 
of rotation and magnetic fields if one is to understand ei¬ 
ther. The strong flow of material from the hydrogen-rich 
envelope, through burning shell(s), and into the core im¬ 
plies that one cannot treat core and envelope in isolation 
from each other. 

This observation leads us to discard models of the ro¬ 
tation which neglect the influence of large-scale poloidal 
magnetic fields that were deposited by previous convec¬ 
tive activity. For example, the Modules for Expe ri ments 
in Stellar Astrophysics (MESA) code (iPaxton et al.l 
1201 111 implements a heuristic prescription for angular mo¬ 
mentum redistribution in the radiative parts of a star, 
which starts with the current-driven (Tayler) instability 
of a predominantly toroidal magnetic field, leading to a 
poloidal component that is then wound back by differen¬ 
tial rotation into the toroidal direction (lSpruitll2002i ). 

The nearly solid rotation of the Sun’s core has been 
used to i nfer the p resenc e of a helical magnetic field 
ijCoueh fe McIntyre! 1199811 . The introduction of a per¬ 
sistent twist to th e fie ld not only facilitates hydromag- 
netic stability (iBraithwaite fc Spruitl2004l ). but also sup¬ 
presses differential rotation across poloidal magnetic flux 
surfaces. The minimal poloidal field that will trans¬ 
port angular momentum rapidly enough to compensate 
changes in the core mass profile on the MS is only ~ ICC 6 
G; but a field ~ 10 2 -10 J times stronger is required to 
compensate a latitude-depenent convective stress that is 
partly transmitted through a tachocline layer into the 
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core (Paper II). The ascent of the early red giant branch 
is completed over a Gyr, not a great deal shorter than 
the present age of the Sun, meaning that a similarly weak 
poloidal magnetic flux would e nforce nearly solid rotatio n 
in the core of a subgiant (e.g. iMaeder fc Mevnetll2f)14l h 

In a 1-2 Mq star, the source of this field differs between 
the early and later stages of post-MS expansion. The en¬ 
velope expands in mass during the first dredge-up phase, 
so that any remnant field in the outer core would have 
been deposited by a vigorously convective and rapidly 
rotating envelope during the pre-MS contraction, when 
magnetic fields several orders of magnitude stronger were 
present. The minimal poloidal magnetic field that pushes 
the core toward solid rotation is well below the threshold 
of detectability if transported to the surface of the star. 
We note that our assumption of solid rotation in radiative 
material depends on the relative weakness of additional 
processes that might source differential rotation, such as 
the damping of gravity waves that are excited in adjacent 
convective zones. 

We are led to a simpler model for the angular mo¬ 
mentum profile in a post-MS star: nearly solid rotation 
in the core, as enforced by large-scale Maxwell stresses, 
combined with inhomogeneous rotation in the convective 
envelope. The details of a dynamo operating around the 
core-envelope boundary, the ensuing flow of magnetic he- 
licity into the core, and the conditions under which the 
dynamo shuts off and core and envelope decouple, are 
addressed in Paper II. Here our focus is on the angular 
velocity profile of the envelope, especially as it is influ¬ 
enced by the Coriolis force. 

Convective Reynolds stresses appear to transport an¬ 
gular momentum inward in accretion flows (IBalbus et al.1 
1996). This effect is demo nstrated in stars with ex tended 
convective envelopes by iBrun fc Palaciosl (|2009H using 
three-dimensional anelastic simulations. In slowly rotat¬ 
ing cases, most relevant to the present discussion, the 
specific angular momentum is found to be independent 
of radius. 

The rotation frequency of the solar convective envelope 
is, by contrast, roughly independent of radius - a fact 
that appears to have influenced treatments of rotation in 
the post-MS phase. Although the solar envelope also cov¬ 
ers many scale heights, it has a modest aspect ratio (0.7:1 
in radius), implying that thermal boundary conditions 
have a stronger influence on the rotation profile than i n 
a star with a deep envelope (e.g. iBrowning et al.1 120061 1. 
The Sun is not especially rapidly rotating by the stan¬ 
dards of late-type MS stars, although the Coriolis force 
still has a significant effect on its rotation profile. The 
rotation period P rot = 2ir/Q is about 3 times the local 
convective time r con = ip/v con (where v con is the speed 
of a convective eddy and £p the pressure scale height), 
corresponding to a Coriolis parameter Co = flr con ~ 2. 
Dynamo act ion may also influence the rotation profile 
(IBrun et al.ll200dll . 

The most important decision to be made here is how 
to handle the backreaction of the Coriolis force on the 
inward pumping of angular momentum. We investigate 
two simple prescriptions: i) a balance between the Cori¬ 
olis force and non-radial pressure gradient forces that 
result from the entropy difference between convective 
downflows and upflows; and ii) an upper limit Co ~ 1 in 
the inner envelope, representing an approximate balance 


between the Coriolis and inertial forces. Both of these 
prescriptions are compared against asteroseismic fits to 
Kepler data for subgiants and post-core He flash stars. 
We find that prescription i) provides a very promising 
quantitative fit. 

Measurements of surface rotation in subgiants pro¬ 
vide independent constraints on the internal an gular mo¬ 
mentu m distribution in post-MS stars. iSchriiver fe Polsl 
(11993T) have used stellar models to argue for an inconsis¬ 
tency between these surface measurements and the evolu¬ 
tion of ~ 1.5 Mq subgiants that maintain solid rotation. 
Adding a more rapidly rotating core, while maintaining 
solid rotation in the envelope, was found to have a neg- 
ligble impact on the result. The discrepency between 
models and data was therefore interpreted in terms of an 
additional magnetic wind torque acting on a star as it 
develops a convective envelope upon leaving the MS. 

Here we show that the same inward pumping of angu¬ 
lar momentum in the envelope that reproduces the Ke¬ 
pler observations also causes a factor ~ 1/2 reduction 
in surface rotation speed over the same interval, thereby 
removing most (but not necessarily all) of the evidence 
for a magnetic wind torque. The calculations presented 
here, and in Paper II, neglect any rotational effect of 
surface magnetic fields during the post-MS phase. 

We also extend our considerations to the tip of the 
RGB and AGB. Here the Coriolis force has a reduced ef¬ 
fect, due to the greater expansion of the star (> 1 0 times 
large r th an the Kepler samp le analyzed by iBeck et ~ahl 
120121 and iMosser et al.ll2012li . This leads us to consider 
more inwardly peaked rotation profiles than the one sug¬ 
gested by fitting Kepler data. 

A final consideration is binary interaction. This is es¬ 
pecially important for stars which lose most of their na¬ 
tal angular momentum to a magnetized wind on the MS, 
corresponding to an initial mass < 1.3 Mq. After such 
a star ascends to the upper RGB or AGB, the inward 
pumping of angular momentum cannot, by itself, sustain 
Co > 1 in the convective envelope. 

A robust hydromagnetic dynamo is maintained dur¬ 
ing the greatest expansion of the star only if it encoun¬ 
ters an external source of angular momentum. We con¬ 
sider the interaction of a converting giant with a com¬ 
panion planet, calculating the combined orbital evolution 
of the planet and the change in the angular momentum 
of the star. An interesting detail of this evolution is 
the excitation of orbital eccentricity by the gravitational 
quadrupole associated with large-scale convective eddies 
in the giant. In the case of a planetary companion, this 
effect is much more pr onounced than when the compan¬ 
ion is a neutron star (lPliinnevlll992f) . and we show that 
it extends the range of orbital periods within which the 
planet gives up its orbital angular momentum to the star. 

1.1. Plan of the Paper 

We investigate the rotation profile of a deep convective 
envelope in Section [2l with a focus on the backreaction of 
the Coriolis force on angular momentum pumping. Var¬ 
ious angular velocity profiles are compared against the 
Kepler measurements of core rotation for subgiants and 
post-He flash stars in Section [3] The effect of angular 
momentum pumping on the evolution of surface rotation 
of subgiants is examined in Section 13.21 

In Section [H we consider the rotational evolution of 
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model stars of initial mass 1 and 5 Mg, including the ef¬ 
fects of convective pumping, angular momentum loss to a 
wind, and planetary interaction. Details of the evolution 
of the orbit of a planet around an expanding giant are 
given in Section [5J including drag from the tide raised on 
the star and the excitation of eccentricity by a convective 
quadrupole. 

Some outstanding questions are addressed in the con¬ 
cluding Section [6] The appendix gives further details of 
our calculation of eccentricity growth in the orbit of a 
companion planet. 

2. PUMPING OF ANGULAR MOMENTUM IN DEEP 
CONVECTIVE ENVELOPES: BACKREACTION FROM 
THE CORIOLIS FORCE 

The deep convective envelope of a giant star divides 
into an outer part that is slowly rotating (Co < 1); and 
an inner part where Co > 1 and the rotation profile is 
influenced by the Coriolis force. 

The entire envelope may rotate slowly near maximum 
expansion, especially if the spin angular momentum of 
the star has been strongly depleted by a magnetized wind 
during the MS phase. On the other hand, a star is still 
compact enough during the early stages of post-MS ex¬ 
pansion that flr con 1 throughout much of the convec¬ 
tive envelope. 

We imagine that entropy is approximately conserved 
over downflows and upflows that are extended com¬ 
pared with a local scale height. Some indirect evidence 
for such a configuration has long been obtained from 
one-dimensional models of giant convection, which re- 
qu ire a m ixing length e xceed ing a single scale height 
dSackmann fe Boothrovdlll99lD . 

Such a suppression of radial mixing is suggestive 
of an inwardly peaked rotation profile in the outer, 
slowly rotating envelope, corresponding to the conser¬ 
vation of specific angular momentum j by individual 
convective flows - as in deed is obtained numerically by 
iBrun fc Palacios! (120091) . Rotation profiles intermediate 
between D(r) oc r~ 2 and solid rotation could clearly be 
obtained from limited mixing, but for definiteness we fo¬ 
cus here on the simple case dj/dr —»• 0 where Co < 1. 
Obtaining the angular velocity profile in the inner enve¬ 
lope depends on a more detailed analysis of the Coriolis 
force, which we give below. 

We consider the evolving partition between fast and 
slow rotation in Section l3Tl and how it depends on the 
angular momentum retained (or gained) by the star in 
Section [4] Given a flat distribution of specific angular 
momentum in the outer envelope, an intermediate-mass 
star generally sustains an inner zone where Co > 1, even 
near the tip of the RGB or AGB; but a star which spins 
down on the MS must gain angular momentum from a 
planetary or stellar companion. 

2.1. Extended Upflows and Downflows with 
Quasi- Geostrophic Balance 

When the inner zone with Co > 1 covers a wide range 
of radius, we can look for a power-law scaling of the ro¬ 
tation frequency Q with spherical radius. The stellar en¬ 
velope is nearly adiabatic and spherically stratified. In 
this inner zone of fast rotation (fllr sin 9 u CO n) about 
an axis z (6 = 0), the steady vorticity equation reads 


(e.g. lBalbusll2009D 

. Bn 2 1 dPdS g(r) dS 

rsme ~d7~ c^~d^~de ~ ~c^~de' 

Here P and p are pressure and density, g{r ) = 

GM(r)/r 2 is gravity, and the specific entropy is S = 
— l)] -1 ln(P/p 7 ), approximated here as that of an 
ideal gas with mean atomic weight p and specific heat 
C p ( C v ) at constant pressure (volume). 

The ‘thermal wind’ approximation used in Equation 
(HD holds because \dP/dO\ -C \dP/dlnr\. In the up¬ 
per, slowly rotating part of the envelope, the latitudinal 
pressure gradient is sourced by the entropy difference be¬ 
tween upflows and downflows, and is smaller by a factor 
(ucon/cs) 2 than the radial gradient. Our focus here is on 
a rapidly rotating solution to Equation 0. where the 
angular pressure gradient which sources dS/dd is pro¬ 
portional to n 2 . 

When the gravitating mass in the inner envelope is 
dominated by the stellar core, gravity g(r) oc r~ 2 . From 
Equation m we obtain the ansatz 

n 2 (r,e) = n 2 (R c )(j^j /( cose). ( 2 ) 

Here R c is the radius inside of which Co > 1. 

One subtlety here is that the stellar mass profile is 
not as centrally concentrated during the early subgiant 
expansion. This generally leads to a shallower scaling 
of it with r: taking g{r) oc r - ' 3 , with /3 < 2, we find 
fl(r) oc r _ ( 1+ ^/ 2 . A fit to the gravity profile in the inner 
convective envelope over the range of expansion that is 
probed by the Kepler asteroseismic data gives /3 ~ 1, 
corresponding to f2(r) oc r — . The monopolar scaling 
of g is approached following the first dredge-up, and is 
maintained near the tips of the RGB and AGB. 

We imagine that the rotation profile organizes at Co ^ 
1 into axially symmetric rolls. In other words, the strong 
latitudinal dependence of the Coriolis force is hypothe¬ 
sized to translate into a persistent latitudinal entropy 
gradient. If, furthermore, the entropy is conserved along 
radially extended convective plumes, then we can decom¬ 
pose the entropy profile into low-order spherical harmon¬ 
ics, S(r, 9 , (j)) = So + SS 0 P e (cos9). Substituting this into 
Equation HD along with Equation 0, and normalizing 
/ = 1 at 9 = 0, we get 

/ = 1 (^ = 2); f = ^ cos 2 9 - (£ = 4 ); 

/=-!- (495 cos 4 9 - 234 cos 2 0 + 19) (£ = 6). (3) 

280 

Angular and radial gradients in Q are generally com¬ 
parable in this type of flow (excepting the simplest case 
of a quadrupolar entropy pattern). We note that rota¬ 
tion is fastest at the pole for a given i (corresponding 
to an ‘anti-solar’ rotation profile), but that the opposing 
gradient can be obtained by superposing £ = 2 with any 
higher harmonic. A small anti-rotation is also present 
in a narrow latitudinal band when the entropy follows a 
single spherical harmonic. 

A caveat here involves our neglect of the Lorentz force, 
which may have a considerable influence near the bound¬ 
ary of a convection zone, but does not obviously compete 
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with the right-hand side of equation 0 within a deep 
convective layer. In particular, the magnetorotational 
instability (MRI) is ineffective in the outer part of the 
envelope where Co < 1, because the growth time of a 
MRI mode is longer than the convective time. 


2.2. Rapid Mixing Between Upflows and Downflows 

We now turn to consider the consequences of rapid 
mixing between upflows and downflows for the rotation 
profile in parts of the envelope that reach Co ~ 1. This 
pro vides a useful comparison with the results of Section 
12.11 and demonstrates how the radial angular velocity 
gradient is tied to the presence of radially extended con¬ 
vective flows. We show that the downward advection of 
angular momentum cannot push the flow to Co 1 in 
the inner part of an extended envelope, where r con de¬ 
creases inward. This result has a simple interpretation: 
the competition between the Coriolis force and the iner¬ 
tial force limits Co to a value of the order of unity. 

In the mixing length approximation, the entropy dif¬ 
ference between upflows and downfl ows is related to the 
speed of a convective eddy by (e.g. iBethanMolf 


_5h_ = C p v 2 on 
k B T C p — C v k B T' 


(4) 


Here Sh is the perturbation to specific enthalpy and T is 
temperature. Then the angular entropy gradient can be 
estimated as dS/d6 ~ (r/£p)SS. Because it turns out 
that Co = 0(1), we cannot take Equation 0 literally, 
but can still use it to estimate the magnitude of the shift 
of angular velocity across a scale height: 


S(fl 2 ) 


9(r) VconV 
Ip k B T 


^con 

ip 


2 


(5) 


We look for a rotation profile that diverges inward, so 
that the constant term in D can be neglected. Then we 
find 

n(r) ~ (6) 

As advertised, the Coriolis parameter does not rise much 
above unity. (Faster rotation can, of course, still be ob¬ 
tained by superposing uniform rotation on the flow.) 

For example, in a nearly adiabatic envelope with 
g{r ) oc r^ 2 , the speed of an eddy scales as u CO n oc 
(pr 2 ) -1 / 3 . Then p(r) oc r -3 / 2 implies f2(r) oc r -7 / 6 for 
7 = 5/3, which is shallower than the profile 0. 


2.3. Model f](r) Profiles 

A simple parameterization of the rotation profile which 
incorporates these lessons may be constructed as follows. 
Where the Coriolis force can be neglected, at r > R c , we 
consider a flat distribution of specific angular momen¬ 
tum, 

n{r) = n(Rc)(^j ; r>R c . (7) 

Here is the angular velocity averaged over a spheri¬ 
cal shell. If tire angular momentum of the star is small 
enough that O(r)r con < 1 throughout the envelope, then 
we identify R c with the base of the envelope, 

Rc — 77b env j (f^T CO n ^ Rb env )■ (8) 


When the star has more angular momentum, we main¬ 
tain the profile 0 in the outer part of the convection 
zone, and consider 

(I) Q(r) = Q(R c )(^j 

COtrans / T \ 

— ( E? \ ( E> ) ’ 4?benv < r < Rc 

r C on{Rc) \RcJ 

(9) 

in the inner envelope. Here Co tr ans = 0(1). As the an¬ 
gular momentum increases R c moves outward, and may 
reach the outer radius of the star, in which case Co is 
larger than unity everywhere in the envelope. 

Following the above discussion, we will explore power- 
law indices 1 < a < 3/2, with the preferred value de¬ 
pending on the radial scaling of gravity. Such a rotation 
profile corresponds to transport of energy by radially ex¬ 
tended, adiabatic plumes. 

A second profile, based on mixing length theory, is ob¬ 
tained by setting Co to a constant everywhere inside the 
radius R c , 

(II) fi( r ) = + S2 0 ; R henv <r<R c . (10) 

The most important difference with profile I is that Co 
remains pegged at Cotrans ~ 1 in the absence of the con¬ 
stant term. Setting flo = 0 strictly limits the angular 
momentum which the envelope may contain, a bound 
which is easily violated during core helium burning. We 
also require that R c is bounded above by the radius at 
which T con reaches a maximum; this prevents the forma¬ 
tion of a layer below the photosphere in which f l increases 
outward. See Figure [I] 

2.4. Latitudinal Differential Rotation 

The latitudinal shift in rotation frequency (e.g. be¬ 
tween pole and equator) has an important influence on 
the growth of magnetic helicity in the growing core of 
an RGB or AGB star (Paper II). In rotation model 
I, where Co > 1 is maintained by radially extended 
convective flows, we deduce a strong latitudinal shift: 
ACl/fl = (S2 eq — Upoie )/Ll ~ 1. A strong shift is also 
maintained in the mixing length approximation (rota¬ 
tion model II), as long as the envelope rotates rapidly 
enough to maintain Co ~ 1 at the base of the enve¬ 
lo pe. It is seen in the mo re rapidly rotating simulations 
of iBrun fe Palaciosl (|2009f ). 

We note that weaker latitudinal differential rotation is 
expected in a compact convective envelope with u con <C 
Qr sin 9, as is encountered in the Sun, Af l/Cl oc Co~ 2 . 

2.5. Matching of Rotation between 
Core and Convective Envelope 

The rotation of core and envelope are easily coupled to 
each other when the inner envelope retains enough an¬ 
gular momentum to sustain a large-scale hydromagnetic 
dynamo. The details of this coupling depend on whether 
the core grows, or recedes, as measured in the radial 
mass coordinate. When the core is growing, the instan¬ 
taneous state of the dynamo is the main consideration, 
because magnetic fields that are amplified near the core¬ 
envelope boundary are advected downward. When the 
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Fig. 1.— Examples of rotation models I (top panel) and II (bot¬ 
tom panel), both taken from the 5 Mq model near the tip of the 
AGB. Each panel shows several rotation profiles corresponding to 
different total stellar angular momenta. The magenta curve shows, 
for reference, the profile of r con (r). Black solid lines represent 
segments where f2(r) oc r ~ 2 and red dotted lines f2(r) oc r~ 1 . 
Top panel: as J* is increased, the transition radius R c , for which 
Co (R c ) = 1, moves outward. Bottom panel: profiles from bot¬ 
tom to top represent models with increasing J*, but now when R c 
reaches the radius of maximum r co n, we beg in to add a uniform 
H component (the last term in Equation (1101) h This prevents the 
formation of an outer zone of positive dQ/dr. Black curve: con¬ 
stant j , reaching Co = 1 at the base of the envelope; green curve: 
fi(r) = v C on/^p', blue curves: constant Q component that grows 
with increasing J*. Kepler asteroseismic data favor model I over 
model II. 


core is receding in mass (as it does when a star first as- 
cends the giant branch during the first dredge-up phase), 
the coupling also depends on a prior dynamo process that 
implanted a magnetic field in the core material. 

A key consideration is the magnitude of the drift speed 
v r — dRhenv/dt of stellar material with respect to the 
core-envelope boundary (averaged if necessary over ther¬ 


mal pulsations, which are present on the terminal AGB). 
Solid rotation is enforced in the outer core during down¬ 
ward drift if the radial Alfven speed is larger than this 
drift speed. This requires that the large-scale poloidal 
magnetic field that is generated in the inner envelope be 
strong enough that 

B? > {'v r - dR henv /dt.) 2 

47rp(i?benv)^ on ~ V* on 

We show in Figure [2] both the drift speed relative to the 
core-envelope boundary, and the convective speed in the 
inner envelope (as evaluated using mixing-length theory), 
during the post-MS evolution of a 1.5 Mq model star. 
The same quantities are also shown for a 5 Mq star near 
the tip of the AGB. (Further details of these models are 
described in Sections [3] and [4]) 

One sees that during the early subgiant expansion of 
the 1.5 Mq model, the remnant poloidal magnetic field 
in the outer core need only retain an energy density that 
exceeds ~ 10 -16 — 10 -13 of the instantaneous convec¬ 
tive energy density in order to compensate changes in 
rotation induced by the evolving mass profile. Magnetic 
fields ~ 10 2 -10 3 times stronger are needed to compensate 
a latitude-dependent convective stress (Paper II). 

The envelope decreases more rapidly in mass near the 
tip of the RGB and AGB, and the magnetic energy 
density that is self-consistently generated during these 
phases must rise to about ~ 10~ 8 —10~ 6 of the convective 
energy density. We conclude that magnetic decoupling 
of core from envelope is most easily achieved during the 
brief phase where the envelope is expelled. 

In Section [5] we concentrate on the comparison with 
Kepler asteroseismic data, which only covers post-MS 
stars more compact than ~ 1O1? 0 . These subgiants 
and core He burning stars evolve on a relatively long 
timescale (~ 10 2 times the duration of the thermally 
pulsating AGB phase). The results shown in Figure [2] 
strongly suggest that the rotation of core and envelope 
remain well coupled in these stars. 

We therefore take the core to rotate as a solid with the 
mean rotation at the base of the envelope, 


fin 


f^benv — f^(Rbenv)- 


( 12 ) 


An explanation for the rotation behavior of the Kepler 
sample must be found in the redistribution of angular 
momentum within an extended convective envelope. 


3. COMPARISON OF 1.5 M© MODEL STAR 
WITH KEPLER DATA 

We begin by calibrating the rotation models described 
by Equations © and ©-dD against the core rotation 
periods of 1 — 2 M © subgiants and p ost- core He flash 
stars a s measured bv lBeck et all (120121 ) and lMosser et al.l 
©Hi). F or ease of comparison , we use the same stellar 
model as iGantiello et all (|2014l ). namely a 1.5 M 0 star 
of solar metallicity with zero-age main sequence (ZAMS) 
equatorial rotation speed 50 km s -1 . 

Our calcul ation also employs the ID stellar evolution 
code MESA (jPaxton et alil20lTL version 5527), but with 
the built-in prescriptions for angular momentum trans¬ 
port turned off. Instead, the rotation profile is ‘patched’ 
onto the sequence of stellar models. 

The total spin angular momentum of the star is evolved 
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Fig. 2.— Flow speed of stellar material with respect to the core¬ 
envelope boundary, compared with convective speed in the inner 
envelope. Upper panel: 1.5 Mq model from the early subgiant 
expansion to the thermally pulsating AGB. Bottom panel: 5 Mq 
model at the tip of the AGB. The logarithmic abscissa in the upper 
panel is referenced to the time £tagb of envelope ejection. The 
time coordinate in the lower panel has been shifted to focus on the 
final few hundred years. 
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self-consistently in response to mass loss, using the pre¬ 
scription described in Section r4.ll Magnetic wind torques 
are neglected, for the reasons described in Section 13.21 
The net effect of mass loss is found to be modest up to 
the onset of core He burning: a factor ~ 2 reduction in 
angular momentum during the RGB phase. 

We observe from Figure[3]that the inner envelope must 
rotate rapidly (Cobenv > 1), even in the case of solid ro¬ 
tation, over the range of expansion probed by the Kepler 
data. The rotation period Pbenv = 27r/flbenv at the base 
of the envelope would actually decrease with increasing 
expansion if the entire envelope could sustain a profile 
fl(r) oc r -2 . 

This effect is still present in the less peaked rotation 



t (years) 


Fig. 3.— Rotation period at the base of the convective envelope, 
in a Mzams = 1.5 Mq model star, as a function of time during 
the RGB and early core He burning phases. Curves correspond 
to different rotation profiles. Black line: rotation period Pbenv 
corresponding to Cobenv = 1 (with r con evaluated in the mixing- 
length approximation). Green dashed line: solid rotation. Red 
dotted line: uniform specific angular momentum (even where Co > 
1). Blue long-dashed line: rotation model I given by Equations 0, 
(lOb with a = 1. Magenta dashed-dotted line: rotation model I with 
a = 3/2. Vertical black lines show the age range probed by the 
Kepler subgiant measurements, with corresponding stellar radius. 


profile I, as given by Equations o, m with a = 1 (or 
3/2). Now Pbenv sits close to the observed range of 10—20 
days. This period is identified here with the core rotation 
period. 

Another important feature is the large jump in rota¬ 
tion period between the red giant and core He burning 
phases (Figure [I}. The latter, slower, rotation results 
from a retraction of the convective envelope combined 
with continuing angular momentum pumping and core- 
envelo pe syn chronization. (It was previously ascribed 
by iCantiello et al.l (|2014| ) to the conservation of angu¬ 
lar momentum by an expanding core that has already 
decoupled rotationally from the inner envelope.) 

A closer look at the effect of the inner rotation index a, 
and the initial rotation speed of the star, is provided by 
Figures 0] and [5] Here one sees that an index a around 
1 — 1.1 provides a close agreement with the Kepler mea¬ 
surements for stars more massive than ~ 1.3M 0 . There 
is some degeneracy between the effects of changing a and 
initial rotation speed. It is interesting to note that this 
profile agrees well with the analytic scaling derived for 
rotation model I in Section O in the case (appropri¬ 
ate to the expansion phase probed by the Kepler data) 
where g(r) oc r _1 in the inner envelope. Recall that this 
rotation model represents radially extended upfiows and 
downflows in a geometrically deep, adiabatic envelope. 

A different conclusion emerges for stars close to so¬ 
lar mass. We see in the bottom panel of Figure []] that 
the core of a 1 Mq star without a close planetary com¬ 
panion is predicted to spin ~ 5 — 10 times more slowly 
during subgiant expansion than is seen in our 1.5 Mq 
model. On the other hand, the core rotation during the 
He burning phase is only ~ 3 times slower, and indeed 
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Fig. 4.— Rotation period at the base of the convective envelope 
during the early giant expansion and the core He burning phase 
of a 1.5 Mq star (upper panel) and a 1 Mq star (lower panel). 
Rotation model I given by Equations 0, {9} with Cotrans = 1 
and varying inner index a (here fl(r) oc r ~“ inside the radius R c 
where Co = Cotrans). The 1.5 Mq star is given a ZAMS equatorial 
rotation speed 50 km s 1 , while the 1 Mq star is given the solar 
rotational angular momentum . Diag onal black line: fit to Kepler 
subgiant data from ICantiello et ali 62014). The radius at which 
majority of helium clump stars are found is marked in the upper 
panel. 

remains within the range measured by Kepler for clump 
stars. Significant spindown by magnetic wind torq ues 
during the subgiant phase (Ivan Saders fc Pinsonneaultl 
12013 ) would have similar effects. 

3.1. Transition from Co > 1 to Co < 1 throughout 
Bulk of Convective Envelope 

The Kepler sample of subgiants and core He burning 
stars probes a different rotational regime than is expe¬ 
rienced by most giants near their maximum expansion. 
The Kepler stars are compact enough that we find Co > 1 
throughout most of the convective envelope, excepting a 



Fig. 5. — Same as upper panel of Figure [T] but now for various 
initial rotation speeds of the 1.5 Mq star, and two selected scalings 
for S2(r) in the rapidly rotating inner zone. 



Fig. 6.— Profile of Coriolis parameter with radius during two 
phases of the evolution of our 1.5 Mq model with initial v TO t = 50 
km s -1 . Black curve: subgiant (stellar radius 6 Rq) with inner ro¬ 
tation index a = 1 in Equation foi l. Colored curves: near the tip of 
the AGB (stellar radius 100 Rq) with uniform specific angular mo¬ 
mentum in outer envelope, and inner rotation profile corresponding 
to either a = 1 or 1.5, or to Co saturating at unity. 


relatively thin layer below the surface. Indeed Co reaches 
a considerable value ~ 10 — 30 at the base of the envelope 
(Figure [6]). The measured core rotation rates probe the 
angular velocity profile that is sustained in the envelope 
in a regime of relatively rapid rotation. 

Figure |T] shows the evolving partition between an outer 
zone with uniform specific angular momentum, and an 
inner zone with angular velocity parameterized by Equa¬ 
tion ©• These curves correspond to a star with spin 
angular momentum close to its birth value, assuming an 
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Fig. 7. — Profile of convective envelope during RGB expansion 
(top panel) and AGB expansion (bottom panel) in rotation model 
I. As the aspect ratio of stellar radius to the base of the envelope 
grows, the domain within which Co > 1 shrinks. The radius R c 
marks the transition point between an outer zone with uniform 
specific angular momentum, and an inner zone with £l(r) oc r ~ a . 
This transition depends only weakly on the inner index a , because 
most of the angular momentum is stored in the outer envelope. 


equatorial surface rotation speed of 50 km s _1 . Out to a 
stellar radius i?.* ~ 10 -Rq, the envelope remains compact 
enough that log(i?*/i? c ) < log(i? c /i?benv)- Beyond this 
size, the inner rotation rate is controlled by the rotation 
profile at Co < 1. 

The spin and magnetization of the WD remnant, as 
explored in Paper II, are strongly influenced by the 17 (r) 
profile of a slowly rotating convective envelope near the 
tip of the RGB and AGB. Then the gravity profile steep¬ 
ens in the inner envelope, and we maintain consideration 
of rotation profiles at Co > 1 that are somewhat steeper 
than suggested by the fit to Kepler data: a ~ 3/2 as 
opposed to a ~ 1 in Equation ([9|). 


M zams =1.5M s base RGB 



Fig. 8. — G reen crosses: observed s urface rotation speeds of class 
IV stars from lSchriiver &; Pols! <119931) . multiplied by 4/7T to account 
for unknown inclination, as a function of effective temperature. We 
show for comparison our model of a 1.5 Mq star on the early red 
giant branch, with initial rotation speed 50 km s —1 . Curves cor¬ 
respond to different profiles of rotation in the envelope, combined 
with solid rotation in the core. Solid rotation is enforced in the 
envelope when the aspect ratio Rbenv/-^* between the base of the 
envelope and stellar photosphere is larger than 0.9 or 0.5. The 
a = 1,1.5 curves correspond to rotation model I. In the case a = 1 
(which best reproduces the asteroseismic fits to Kepler data) there 
is a factor ~ 1/2 reduction in surface rotation compared with a 
solid rotator at log(T e ff) <3.8. 


3.2. Surface Rotation of Subgiants: 

Angular Momentum Pumping Versus Magnetized Winds 

The expansion of a star following the MS phase 
is associated with a drop in surface rotation. The 
strength of this drop is sensitive to the distribution 
of angular momentum within the star, as well as to 
external torques associate d with magnetic activi ty in 
t he deepening en velo pe. S chriiver fc Pols (ll993t l and 
Ivan Saders &; Pinsonneaultl (12013H have presented evi¬ 
dence for an inconsistency between solid rotation in the 
envelope and measurements in the surface rotation of 
1.5—2 Mq stars: the measured rotation is a factor ~ 2—3 
slower than is expected for constant stellar angular mo¬ 
mentum. 

A similar effect would be imposed by inward p urnp- 
ing of angular momentum. iSchriiver fc Polsl (I1993D show 
that a fast-spinning core has a negligible effect on the re¬ 
sult if the convective envelope maintains solid rotation, 
but did not consider envelopes with inward-peaked ro¬ 
tation profiles. The key point here is that the radiative 
core becomes very compact early in the expansion onto 
the red giant branch, in the sense that its moment of 
inertia is a very small fraction of Mb env R^ env . There¬ 
fore inhomogeneous rotation in the envelope has a much 
larger effect on surface rotation than does a fast-spinning 
core. 

We show in Figure [8] the evolution of the surface ro¬ 
tation in our 1.5 Mq model, for the same set of inho¬ 
mogeneous rotation profiles that we have previously ex¬ 
plored. The envelope at high T e g is still relatively thin. 
We therefore impose <9I7/i dr = 0, to mimic the solar rota- 
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tion profile, when the aspect ratio i?benv/R* rises above 
a critical value 0.5 or 0.9. The profile that is favored 
by the Kepler core rotation data, Equations © and © 
with inner index a = 1 , shows a factor ~ 1/2 reduction in 
surface rotation compared with the solid rotator. This 
removes a significant p art o f the discrepency found by 
iSchriiver fe Polsl (1199311 and Ivan Saders fe Pinsonneauld 
(12013 1. but perhaps not all. 

4. EVOLVING ROTATION PROFILE WITH MASS LOSS 
AND INTERACTION WITH A PLANET 

We now consider the rotational evolution of a star as it 
approaches the tip of the RGB and AGB. Here the prob¬ 
ability of interacting with a companion increases signif¬ 
icantly. The Coriolis parameter also drops markedly at 
the base of the convective envelope. We maintain our 
focus on a two-layered rotation model. 

Two test cases are analyzed in detail: a star with initial 
mass Mzams = 1 M 0 , which experiences strong spin- 
down on the MS and interacts with a planet; and an 
intermediate-mass star (Mzams = 5 M 0 ) which remains 
rapidly rotating at the end of the MS and leaves behind 
a massive WD. This second stellar model is also allowed 
to interact with a companion. MESA is used to evolve 
these stellar models (both of solar metallicity) from the 
ZAMS to the post-AGB phase. 

In the case of the solar-mass star, the interaction with 
a planet was followed from the base of the RGB, when 
the star has expanded to a radius R* ~ 1O.9R 0 (age 
tio .9 ~ 12.37 Gyr). At that point, the spin angular mo¬ 
mentum was set equal to that of the present Sun. Stellar 
quantities at times intermediate between the outputed 
MESA models were obtained by linear interpolation. De¬ 
tails of the orbital integration and the interaction be¬ 
tween planet and star are presented in Section [5] 

The equatorial rotation speed of the 5 M 0 star is set 
to 50 km s -1 at the ZAMS, correspo nding to the pea k 
in the measured v sin i distribution ([Wolff et al.l 120071 1. 
Both the 1 Mq and 5 M 0 models are taken to rotate 
as a solid body on the MS. Deviations from solid body 
rotation develop after the formation of a deep convective 
envelope, following the approach of Section [5] 

There are instances around the tip of the AGB when 
the outer envelope splits into several convection zones, 
separated by radiative layers. The entropy gradient in 
these radiative zones, although positive, remains small. 
We impose a matching of specific angular momentum 
across them. 


0.1 
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Fig. 9. — Evolution of spin angular momentum of 1 Mq model 
star during RGB, core He burning and final AGB phases. Curves 
represent cases described in the text. Red lines (Jms) : no exter¬ 
nal source of angular momentum. Black lines (./jupiter,^ =iau) : 
absorption of a Jupiter-mass planet orbiting initially at a$ = 1 
AU. Short dashed lines = c): solid rotation throughout core 
and envelope. Other lines correspond to two-layer rotation model, 
with uniform specific angular momentum in outer envelope, and 
varying inner rotation profile. Solid (long-dashed) curves: a = 1 
(a = 1.5) in Equation J9]); long-dashed-dotted curve: profile (1101 ) 
with Co = 1 in inner zone. 
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with r)R = 0.5. Here (~) denotes stellar parameters mea¬ 
sured in solar units. On the AGB we use 

M* Bloecker = 1-93 X 10~ 21 7 Jb ^ Mq yi' _1 (15) 

M* 

from 1 Blocked (il995l) . Mass loss according to the Blocker 
formula is turned on when the central H and He abun¬ 
dances are less than 10 -2 and 10 -4 , respectively. 

The choice of AGB mass loss parameter rin is m oti- 
vated by recent cluster data (e.g. iDobbie et al.ll2009l f. As 
we discuss further in Paper II, we set r/B = 0.05, corre¬ 
sponding to a final WD mass Mwd ~ 0.87M 0 (from the 
5 Mq progenitor). The same wind parameters are used 
for both the 1 M 0 and 5 M 0 models. The solar-type star 
loses ~ O.25M 0 on the RGB and another ~ O.21M 0 on 
the AGB for a final WD mass Mwd ~ O.54M 0 . 


4.1. Loss of Angular Momentum on the RGB and AGB 

Magnetic wind braking is not taken into account in 
our simulations during the subgiant or giant phases, for 
the reasons given in Section [3© Then the ejected mass 
carries the specific angular momentum of the stellar pho¬ 
tosphere, 

■h = ^f2(R*)R*M*. (13) 

The rate of mass l oss o n the RGB is handled using the 
prescription of lReimersl (119751) , 

L R 

M* Reimers = 4 X 10~ 13 7lR *- * 

M* 


4.2. Rotation Rate of the Inner Envelope and Core 

There is a simple relation between spin angular mo¬ 
mentum and core angular velocity Dbenv when the enve¬ 
lope is extended and maintains uniform angular momen¬ 
tum per unit mass, 

j(r) = ^M^j = \ r2 ^ r) = const (16) 

Given a net stellar angular momentum J*, and an effec¬ 
tive moment of inertia in the convective envelope, 

IeS = y j 47rr R henv p(r)dr = — R benv (M* — M\, env ), 

(17) 


Mq yr-\ (14) 










































10 




tip RGB M zams =M q 
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Fig. 10. — Evolution of the Coriolis parameter at the base of 
the convective envelope. Curves correspond to those in Figure [9] 
Upper panel shows the RGB, core helium burning phase, and AGB. 
Lower panel shows an expanded view of the tip of the RGB and 
the beginning of the core He burning phase. 


we obtain 


left + I core. 

Excepting near the point where the envelope has been 
ejected at the end of the AGB, the moment of inertia of 
the core can generally be neglected, I cole <C I e g. 

The angular momentum stored in the core is enhanced 
by a factor ~ {R*/Rbenv) 2 , as compared with the case of 
solid rotation throughout the star. This simple rotation 
model breaks down when it predicts Co > 1 at the base 
of the convective envelope, and we employ the modified 
rotation profiles ©, Cod. 

Considering first the 1 M© model, we compare several 
scenarios. The star is either assumed to evolve in isola¬ 
tion; or to interact with a Jupiter-mass planet that starts 
in a circular orbit of semi-major axis 1 AU. The internal 



t (years) 

Fig. 11.— Effect on the 5M© model star of the engulfment of a 
substellar companion of mass (3,10,30,100)Mj up i te r from an initial 
orbit di = 2 AU. Orbital angular momentum of planet is added to 
star when it reaches a radius R* = 200R© on the AGB, as guided 
by more detailed tidal evolution calculations. The physical con¬ 
sequences depend weakly on the epoch of engulfment, if it occurs 
before the first thermal pulse (at t ~ 1.0717 x 10 8 years), when the 
envelope dynamo begins to contribute to the magnetic field of the 
hydrogen-depleted core. Upper panel: time evolution of transition 
radius R c between outer slowly rotating zone (Co < 1) and inner 
rapidly rotating zone (Co > l). Bottom panel: change in Cobenv 
due to the added angular momentum. Rotation profile corresponds 
to a = 1 in Equation ©■ Although R c increases with increasing 
angular momentum, Cobenv decreases. 


rotation profile is either described by a two-layer model 
(Section [OP during its post-MS evolution (‘j = c’ in the 
outer envelope, with various inner f2(r) profiles); or it is 
assumed to be solid-body from the ZAMS to the terminal 
AGB = c’). 

Figure [9] shows the evolution of J*. We see that solid- 
body rotation leads to a significant loss of angular mo¬ 
mentum through the RGB winds, because little angular 
momentum is stored deep in the star. Introducing an 
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Fig. 12. — Same as bottom panel of Figure eh but now with 
steeper inner rotation profile a = 1.5 in Equation ©• Now Cobenv 
increases monotonically with the added angular momentum. 


outer zone of uniform specific angular momentum allows 
only 1/2 of the angular momentum present at the end 
of the MS to be lost, with the remainder retained at the 
onset of core He burning. The engulfment of the Jupiter 
greatly augments J*, but if 11 is uniform then most of this 
additional angular momentum is lost to RGB winds. 

The most important consideration for us is whether 
Cobenv reaches the threshold for dynamo activity. The 
time dependence of Cobenv is shown in Figure [TOl Only 
in the case of planet absorption with ‘j = c’ in the outer 
envelope does the star comfortably pass this threshold. 
Enough angular momentum is absorbed from a Jupiter- 
rnass planet that almost any profile H(r) oc r~ a with a > 
0 will allow the star to attain Cobenv = 1 on the RGB. 
Stronger angular momentum pumping is required near 
the tip of the AGB to sustain rapid rotation: uniform 
rotation combined with absorption of a Jupiter does not 
suffice. 

An isolated -Mzams = 1 M 0 star without a planetary 
companion does not attain Cobenv ~ 1 in its inner en¬ 
velope near the tip of the RGB or AGB, no matter the 
rotation profile. This dynamo threshold is missed by 
more than two orders of magnitude in the case of solid 
rotation. 

Angular momentum could also be transferred to an 
evolved giant by a tidal interaction with a stellar com¬ 
panion. This happens inevitably during the evolution 
to a common envelope and the formation of pre-CV 
system. For example, the 5 Mq star may easily engulf 
a companion during the AGB phase, given the large 
maximum radius reached by the star (~ 6OOi? 0 ). We 
consider the effect of injecting a companion of mass 
(3,10,30,100)Mj U p it er from an orbit with initial semi- 
nrajor axis cq = 2AU when the star reaches a radius 
i?+ = 2OOi? 0 on the AGB. We also consider a range of 
inner rotation profiles, a = 1 to 1.5 in Equation ([9]) . 

Figure |TT] shows that the added angular momentum 
pushes R c outward, as expected; but when a = 1 the 
relative flatness of the r con (r) profile (r con increasing 


roughly in proportion to r) may actually cause Cobenv to 
decrease. In other words, the engulfment of a more mas¬ 
sive companion may result in slower rotation at the base 
of the convective envelope. Increasing the rotation index 
to a = 1.5 allows Cobenv to grow with increasing injected 
angular momentum (Figure [12]). Our inference that an¬ 
gular momentum pumping in the envelopes of subgiants 
can sustain Cobenv ~ 10 (see Figure 0 suggests that the 
rotation profile in the inner envelope may be significantly 
steeper than a = 1 near the tip of the RGB and AGB. 
This would imply an even stronger magnetization of the 
WD remnant, as we show in Paper II. 


5. ORBITAL EVOLUTION OF A PLANETARY 
COMPANION TO A 1M 0 STAR 

Our calculation of the interaction of a planet with an 
evolving giant star follows the orbit in some detail. We 
take into account the dissipation of the tide raised on 
the star, and gravitational scattering off a fluctuating 
quadrupole moment. The second effect, which has been 
previously co nsidered in the context of massive neutron 
star binaries (lPhinnevlll992r) , has a weaker radial depen¬ 
dence than the tidal torque and therefore modestly in¬ 
creases the range of semi-major axis over which a strong 
interaction between planet and star can take place. 

We evolve the planet’s orbit in response to these ex¬ 
ternal forces, considering a range of masses (Earth, Nep¬ 
tune or Jupiter) along with different initial semi-major 
axes (at = 0.5,0.75,1,1.5, 2 AU) and eccentricities (e.; = 
0,0.1,0.5). The change in orbital angular momentum is 
deposited in the star and tracked self-consistently, includ¬ 
ing the effect of stellar mass loss described by Equation 

m- 

The tidal acceleration is calculated using (|Hutlll98lll 


d 2 r GM p (H, 

w = 


(1 + 3-r)f - (n env - u)rq 
r 


. ( 19 ) 

In the above expression, the position of the planet is fol¬ 
lowed in spherical coordinates, with u its orbital angular 
velocity and M p its mass. We set the rotation frequency 
fl env of the stellar envelope^ to zero. The time lag of the 
tidal bulge is given by 


where 


Ri 

GM*r con 

(20) 

M env Riy s 

(21) 


is the global convective overturn time. Th e Love number 
fc? is obtained (see Eq uations (11) of Zahn 1989] and (17) 
of I Schar lemand ll 98 2j) . 

D 1 / 3 pi 

k 2 = 20.5a 4 ' 3 / x 22 ' 3 L 4 ' 3 p 2 ' 3 lpdx. (22) 

Menv J x benv 

Here a = 2 is the convective mixing length parameter, 
M env = M* — Mbenv is the mass of the convective enve¬ 
lope of the star, x = r/R+, ccbenv = Rbenv/R*, and L r is 
the luminosity at a given radius. 


1 Equation m is derived assuming solid rotation in the enve¬ 
lope. 
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Fig. 13.— Red line: orbital semi-major axis of a Jupiter-mass 
planet interacting with a RGB star and starting from a circular 
orbit with a,i = 1 AU. Black line: stellar radius. Spiral-in of planet 
just below the tip of the RGB deposits angular momentum and 
triggers a dynamo process near the core-envelope boundary (Paper 
ii). 
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Fig. 14.— Expanded view of the orbital parameters of the 
Jupiter-mass planet in the same realization shown in Figure 1 13 1 
Engulfment of the planet is preceded by a stochastic growth in or¬ 
bital eccentricity, due to gravitational scattering off the convective 
quadrupole. 


The acceleration from the convective quadrupole is cal¬ 
culated using the standard expansion of the gravitational 
potential (see [Murray fc Dermottlfl999l ). 


v —<L 

v quad ~ r 


2r 5 . 

xx 
6(Iyz yz 


(J-yy T Izz 2 Ixx)x 


H - (1-XX Izz 2 Iyy)y H - {I XX H - lyy 2 I z z)z 


(23) 
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I xy xy) 


(In practice we evolve the orbit of the planet in Carte¬ 
sian coordinates.) The moment of inertia tensor Rj is 
normalized using the rms convective quadrupole 


Qc 


~ 7.41 x 1CT 3 / 


M Rm Rl 


r 2/3 R 5/3 
^*37-°-* 13 


(M*/M 0 )(M env /M 0 ) 2 / 3 ' 
(24) 

Here / is a factor of order unity which depends on the 
stellar structure. This expression is obtained by integrat¬ 
ing the density perturbation caused by convection inside 
the star over the envelope; see the Appendix [X] for a 
derivation. 

At time intervals separated by r con , the diagonal 
components of the moment of inertia matrix for the 
quadrupole were chosen randomly, with a magnitude uni¬ 
formly distributed in the range [0, Q C on]- The remaining 
components were then obtained by a random rotation of 
the diagonal matrix. Given the random nature of this 
forcing effect, we ran ten realizations for each set of ini¬ 
tial conditions to obtain a statistical sample of results. 

A planet that is exposed to the tidal acceleration m 
and quadrupolar potential (1231) will evolve in semi-major 
axis and eccentricity. It can either collide with the en¬ 
velope of the star or be ejected from the system. Here 
we are interested primarily in initial conditions leading 
to a collision. In such an instance, we assume that drag 
and dynamical friction act instantaneously, resulting in 


the planet’s engulfment and the transfer of its orbital 
angular momentum to the giant. 

5.1. Jupiter-mass Planet with m = 1 AU 

Consider a Jupiter-mass planet starting at a* = 1AU 
and ti = 0. In Figure [13] we show one realization of 
its orbital evolution, comparing semi-major axis a with 
R+. In this case the planet is engulfed shortly before 
the star reaches the tip of the RGB. We treat this as our 
fiducial model for a Mzams = 1M© star interacting with 
a planet. 

The corresponding eccentricity evolution is shown in 
Figure [TT] The stochastic nature of the convective forc¬ 
ing of orbital parameters is evident here. 

5.2. Other Planetary Configurations 

Figure [15] shows the probability of a planet’s engulf¬ 
ment as a function of the initial (specific) orbital angu¬ 
lar momentum. Results for the three planetary masses 
we consider are shown separately. The separate effect of 
shutting off the convective quadrupole (Q co n —> 0) is also 
shown. 

The convective forcing has an especially strong effect 
on the probability of engulfment for Earth-mass planets. 
In the case of a Jupiter-mass planet, there is a ~ 20% 
enhancement in the range of semi-major axes that re¬ 
sults in a merger. These effects arise mainly from the 
excitation of eccentricity and a corresponding reduction 
in periastron. 


6. CONCLUSIONS 

We have considered the inward pumping of angular 
momentum in a deep stellar convective envelope, with a 
goal of defining a workable rotation model for giant stars. 
Independent theoretical and observational lines of argu¬ 
ment point to a profile that differs radically from those 
obtained from stellar evolution codes such as MESA, in 
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Fig. 15. — Probability of engulfment of planets of various masses 
(Earth, Neptune, Jupiter), as a function of the initial specific or¬ 
bital angular momentum. Negative-sloping dash: planet scatters 
gravitationally off convective quadrupole of giant star. Positive- 
sloping dash: Q CO n = 0. Probability of engulfment is enhanced by 
interaction with the quadrupole, especially for lower-mass planets. 


which the implementation of angular momentum trans¬ 
port typically results in solid rotation within convective 
layers. The work presented here and in Paper II also 
demonstrates that much improved dynamo models must 
be implemented in order to capture basic features of stel¬ 
lar rotation, such as the relaxation to solid rotation in 
radiative parts of a star, and the rotational coupling be¬ 
tween a radiative core and the base of an extended con¬ 
vective envelope. 

The rotation of evolving stars cannot be understood 
even in qualitative terms without including the magnetic 
field built up by a hydromagnetic dynamo operating near 
the boundary between convective and radiative zones. 
Some features of this dynamo are special to giant stars. 
First, the high radiative energy flux leads to a strong 
broadening of the tachocline layer situated beneath the 
core-envelope boundary, which may facilitate mixing of 
core material into the envelope. The growth of the core is 
also accompanied by an inward flux of magnetic helicity. 
These effects are investigated in detail in Paper II, where 
we motivate a strong magnetic coupling between core and 
inner envelope when Cobenv > 0.1 — 0.3. 

When this condition is satisfied, a simple picture 
emerges: nearly solid rotation in the radiative parts 
of the star that previously supported a dynamo, com¬ 
bined with an inhomogeneous rotation profile in a deep 
convective envelope, whose slope depends on the local 
Coriolis parameter. For example, solid rotation in a ra¬ 
diative layer adjacent to convective zone is facilitated 
by a tiny helical magnetic field. Differential rotation 
that is sourced by a changing stellar mass profile dur¬ 
ing subgiant expansion can be significantly reduced if 
the poloidal magnetic pressure exceeds ~ 10~ 16 -10~ 15 
of the convective pressure, increasing to ~ 10~ 8 during 
the thermally pulsating AGB. Magnetic fields ~ 10 2 -10 3 
times stronger are required to compensate a latitude- 


dependent convective stress that is partly transmitted 
through a tachocline layer (Paper II). 

A magnetic field exceeding this minimal value is plau¬ 
sibly present in the radiative core of a subgiant: erasure 
to a lower level would require a very efficient mechanism 
for eliminating magnetic helicity. Inhomogeneous rota¬ 
tion could, in principle, be maintained by a supplemen¬ 
tary mechanism of angular momentum transport such as 
gravity wave damping. 

Our results can be summarized as follows. 

1. Strong inward pumping of angular momentum 
is well motivated by me asurements of core rotation in 
stars of radius < !Oi? 0 dBeck et ahll2012t iMosser e t al 


2012), and by numerical simulation (Brim & Palacios 


200911 . A rotational decou pling between core a nd en¬ 
velope (as envisaged e.g. bv lCantiello et al.ll2014f l is dif¬ 
ficult to achieve over the long duration of subgiant ex¬ 
pansion, depending on an unrealistically small poloidal 
magnetic field. We have found that a profile D(r) oc r -1 
within parts of the envelope that are very rapidly ro¬ 
tating (Co > 1), combined with solid rotation in the 
radiative core, provides an excellent fit to the data for 
stars more massive than ~ 1.3 Mq. This corresponds to 
the profile that is derived from quasi-geostrophic balance 
in an adiabatic envelope with heat transported by radi¬ 
ally extended convective upflows and downflows that do 
not mix significantly with each other over a scale height. 
Given a gravity profile g(r ) oc with /3 ~ 1, we deduce 
f l(r) oc r~ 


More focused measurements of individual stars can in 
principle be used to test this rotation model in partic¬ 
ular, to constrain the range of fI within the convective 
envelope - and also to test for the existence of strong dif¬ 
ferential rotation across the cor e-envelope boundary, as 
in the model of lCantiello et all (|2014f ). Efforts so far to 


measure the full radial profile D(r) have considered sep¬ 
arate zones of solid rotation in core and envelope, as well 
as smoother profiles; the fits are mainly sensitive to t he 
rotation rate in the core (|Deheuvels et al.ll2012t l20l4l) . 

Solar-mass stars which spin down significantly on the 
MS, and do not absorb a planetary compansion, are pre¬ 
dicted to retain much slower core rotation during sub¬ 
giant expansion, but to rotate within the range of rates 
inferred from Kepler observations of core He burning 
stars. Although the Kepler subgiant sample is mainly 
composed of stars more massive than 1.3 Mq, this pro¬ 
vides a sharp test of the rotation model presented here. 
On the other hand, we note that convergence in the core 
rotation of stars below and above the 1.3M 0 threshold 
during subgiant expansion would be remarkable in al¬ 
most any rotation model, given the significant difference 
in rotation rate at the end of the MS. 

2. Independent constraints on internal rotation in sub¬ 
giants come from measurements of the slowdown in sur¬ 
face rotation as stars of mass > 1.3M 0 first ascend the 
giant branch. A good p art of the discrepency, found by 
iSchriiver fe Polsl (119931 ), with stellar models which as¬ 
sume solid rotation, is removed by introducing the same 
envelope rotation profile that is implied by the Kepler as- 
teroseismic data. A separate rapidly rotating core, com¬ 
bined with solid rotation in the envelope, has a negligible 
effect on surface rotation because / cor e is a very small 
fraction of (and of the effective moment of 
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inertia of the envelope). 

3. Stars near the tip of the RGB and AGB (radii 
> 10 2 Rq ) probe a different regime of slow rotation that 
is not directly constrained by the Kepler observations. 
Here we have considered steeper rotation profiles, j ~ 
const, in the outer parts of the envelope where Co < 1. 
This stronger level of angular momentum pumping al¬ 
lows rapid rotation to be sustained deep in the envelope 
even during the final stages of envelope ejection. It also 
limits the loss of angular momentum to winds on the 
RGB and AGB. We find that the conditions for a hydro- 
magnetic dynamo are favorable in a star which rotates 
rapidly at the end of the MS (e.g. initial mass >1.3 M 0 ); 
or in a solar-mass star which absorbs a Jupiter (or even 
Neptune) at some point during its post-MS expansion. 

The gravity profile steepens significantly in the inner 
envelope following the first dredge-up phase, and ap¬ 
proaches g{r) oc r~ 2 near the tips of the RGB and AGB. 
If the convective upflows and downflows remain radially 
extended, then our scaling solution to the vorticity equa¬ 
tion implies that the inner angular velocity profile steep¬ 
ens to Q(r) oc r -3 / 2 . 

4. The mass of an individual convective cell in a gi¬ 
ant envelope greatly exceeds the mass of Jupiter. This 
means than the orbital eccentricity of any planet will 
be significantly excited by the associated gravitational 
quadrupole. We have calculated the orbital evolution 
of a suit of planets of various masses (Earth, Neptune, 
Jupiter) in response to this gravitational stirring as well 
as tidal drag. The gravitational stirring has a significant 
effect on the probability of absorption of an Earth-mass 
planet starting at a semi-major axis 1 AU. 


5. The case of a solar-mass star which spins down 
on the MS and then does not regain angular momentum 
from a companion is more complicated. The rotation 
model described here can be applied while the star re¬ 
mains more compact than i?* ~ 10 12 cm; but slow ro¬ 
tation is enountered throughout the convective envelope 
near the tip of the giant branches. In these most ex¬ 
tended phases of expansion, the coupling between core 
and envelope will depend on additional (and perhaps 
largely hydrodynamic) mixing mechanisms, whose exact 
nature is still be ing debated (IMaeder fc Mevnell 120001 : 
iFuller et al.ll2014ll . 

More extended simulations of deep convection, in¬ 
cluding radiative transport and feedback from magnetic 
fields, are required to confirm the dependence of rotation 
profile on Coriolis parameter. A more centrally peaked 
profile is plausibly achieved near maximum expansion 
on the giant branches, where the effects of rotation are 
weaker, but whether this profile really approaches uni¬ 
form specific angular momentum needs further testing. 
We emphasize the degeneracy between the combined ef¬ 
fects of flattening the slope of the rotation profile and in¬ 
creasing the spindown on the MS (both of which reduce 
the rotation rate in the inner envelope), and raising the 
angular momentum that is received from a sub-stellar 
companion during post-MS expansion. Although mag¬ 
netic field amplification in the inner envelope becomes 
more difficult as the rotation profile flattens, even a plan¬ 
etary companion can have a significant compensating ef¬ 
fect. 

This work was supported by NSERC. We thank the 
anonymous referee for detailed comments on the presen¬ 
tation of our results. 


APPENDIX 

NORMALIZATION OF CONVECTIVE QUADRUPOLE DURING GIANT EVOLUTION 

The internal density, temperature and luminosity profiles of a model giant star evolve with time in a complicated 
way. Here we describe a simple prescription for the net gravitational quadrupole that is induced by convective motions 
in the giant envelope, which can be expressed in terms of integral quanties such as stellar luminosity L*, radius i?* 
and core and envelope masses M core , M env . 

We estimate the magnitude of the quadrupole from 

Qcon ^ (SMr 2 ) (Al) 

where Weils is the number of convective cells in the envelope of the star, and SM is the mass perturbation induced 
by convection in a cell, and the average is over the volume of the envelope. We assume that the convective cells (of 
angular size AQ) penetrate the whole star, so that 

Weils ~ (6Mr 2 ) ~ AS! J r A Spdr. (A2) 

The density perturbation Sp ~ pn 2 on /c 2 , where we estimate pi> 3 on ~ L*/47rr 2 . Therefore 

/ T \ 2 / 3 r R * n 1 / 3 

(SMr 2 ) ~ AQ J r 8 ^ 3 — 2 _c ^ r > (A3) 


where we take AQ = 1. 

The assumption of a point gravitational potential, in combination with an adiabatic equation of state, gives a good 
fit to the actual density profile of a giant. The temperature profile is more sensitive to the mass distribution in the 
star. Hence 


p{r) 


4M env //?.* A 3/2 
n 2 R 3 I, r J ’ 
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where M(< r) is the enclosed mass, which is proportional to 

Mcore , 16M env 


g{x, M core /M*) = 


M* 


ttM, 


* Jo 


x' 2 ( — — 1 ) dx . 


3/2 


The quadrupole m works out to 

Q con 3 AQ 1 / 2 Ll /3 Rl /3 


M env R J 4V37T5/6G- 


m 2 1 3 m+ 


/; / = 


„8/3 


'o 5 (x,M core /M*) 




- 1/2 


dx. 


(A5) 


(A6) 
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